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ABSTRACT 

Context. Filamentary structures are common in molecular clouds. Explaining how they fragment to dense cores is a missing step in 
understanding their role in star formation. 

Aims. We perform a case study of whether low-mass filaments are close-to hydrostatic prior to their fragmentation, and whether their 
fragmentation agrees with gravitational fragmentation models. For this, we study the ~6.5 pc long Musca molecular cloud that is an 
ideal candidate for a filament at an early stage of fragmentation. 

Methods. We employ dust extinction mapping in conjunction with near-infrared JHK^ band data from the CTIO/NEWFIRM instru¬ 
ment, and 870 pm dust continuum emission data from the APEX/LABOCA instrument, to estimate column densities in Musca. We 
use the data to identify fragments from the cloud and to determine the radial density distribution of its filamentary part. We compare 
the cloud’s morphology with *^CO and C'*0 line emission observed with the APEX/SHeEl instrument. 

Results. The Musca cloud is pronouncedly fragmented at its ends, but harbours a remarkably well-defined, ■'-1.6 pc long filament in 
its Center region. The line mass of the filament is 21-31 Mq pc“' and FWHM 0.07 pc. The radial profile of the filament can be fitted 
with a Plummer profile that has the power-index of 2.6 ±11%, flatter than that of an infinite hydrostatic filament. The profile can also 
be fitted with a hydrostatic cylinder truncated by external pressure. These models imply a central density of -“S-lOxlO** cm“^. The 
fragments in the cloud have a mean separation of ~0.4 pc, in agreement with gravitational fragmentation. These properties, together 
with the subsonic and velocity-coherent nature of the cloud, suggest a scenario in which an initially hydrostatic cloud is currently 
gravitationally fragmenting. The fragmentation has started a few tenths of a Myr ago from the ends of the cloud, leaving its centre yet 
relatively non-fragmented, possibly because of gravitational focusing in a finite geometry. 
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1. Introduction 


Filamentary structures have been known to be common mor¬ 
phological fea tures of the i nterstellar medium (ISM) for 
decades (e.g., iBarnardI I1927I: [ Schneider & ElmegreenI I1979I: 


iDobashi et al.l2005l:lAndre et al.i2014h . Observatiorial works ex¬ 
amining the physical properties of filaments have built a pic¬ 
ture in which filamentary clouds fragment into denser clumps 
and/or de nse cores that can further collapse to form new 
stars (e.g., Ballv et al.lll987l: Mizuno et al.l 1995I: Mveri 2009t 


Andr6 et alT 20ldr 2014 ~|.Tackson et al.l l2010t ISchneider et all 

201(i Kainulainen et al.l 20131) . This phenomenology is appar¬ 
ent across a wide range of size- and mass-sca ks, spanning ~0.1- 
100 pc in size and ~1-10^ M, 


(e.g..lPineda et al 


Hacar & Tafalla [201 it Schneider et al.l l2010t iHill et al 


201C; 


2011 


Hacar et al.ll2013tlGoodman et al.l^ldtlRagan et al.ll20i4 ). The 

new, sensitive dust emission observations by the Herschel Space 
Observatory have also manifested the universality of filamentary 
structures in molecular clouds, reinforcing a hypothesis that they 


* This publication is based on data acquired with the Atacama 
Pathfinder Experiment (APEX). APEX is a collaboration between 
the Max-Planck-Institut fur Radioastronomie, the European Southern 
Observatory, and the Onsala Space Observatory (Max-Planck pro¬ 
gramme ID M-085.F-0027). 


may provide a dominan t pathway to star fo rmation in molecular 
clouds (for a review, see lAndre et al.l20T^ . Explaining how fila¬ 
mentary structures evolve and form stars is then a necessary step 
in understanding the exact role of this pathway. 


Analytical studies have established a viable framework for 
the existence of filamentary structures, regardless of their ori¬ 
gin, a s hydrostatic, self-gravitating, thermally pre ssurised cylin¬ 
ders (lOstrikeil 119641 : iFischera & Martiril 12012ah . These cylin- 
ders can fragment in time-scales of ~Mvr when perturbe d 
(llnutsuka & Mivarri^ 1 19921 119971 : IFischera & Martiril l2012all . 
and even small density perturbations in the fila ment can collaps e 
faster than the global collapse takes place in it dPon et al.ll201 ll) . 
This gravitational fragmentation provides one viable model to 
explain the ’’dense cores” observed in filaments. The details of 
the fragmentation process in this model depend most crucially 
on the line mass of the filament at its initial condition, i.e., 
at the equilibrium configuration. The line mass at equilibrium, 
and hence the fragmentation process, can be affected by the 
possible pressure confinem ent by the surrounding environment 
( Fischera & Marti j2012a ). support provided by magnetic fields 
( Fiege&^udritz 2000allb ). accretion of gas onto the fila ment 
( Heitschl 201^ b|), temperat ure profile (iRecchi et al.ll2Q13l) . and 
rotation ( Recchi et al.ll2014l). F or a more comprehensive review, 
we refer to Andre et al. 120141) . 
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One key open question related to filaments is whether i) the 
filaments are close to hydrostatic objects prior to their frag¬ 
mentation, in which case they first have to detach from the 
turbulent regime of the cloud and then fragment, or ii) the 
evolutionary picture of filaments is dynamic and hydrostatic 
conditions are never prevalent in the cloud. This question can 
be addressed by first carefully identifying filaments that are 
at the earliest stages of their evolution, possibly not yet frag¬ 
mented, and by comparing their structure with the predictions 
of the hydrostatic equilibrium models. However in practice, do¬ 
ing this is severely hampered by the fact that non-fmgmented 
filaments that carry the imprints of the initial condition are 
rare. This alone indicates that the phase in which a filament 
is non-fragmented and yet CO-emitting (so that it would be 
identified as a molecular cloud) is relatively short. While sev¬ 
eral recent works have made efforts to estimate the structure 
of filaments in nearby clouds (e.g., Arzoumanian_gLaL 20]Tt 
Fischera & Martini l2012bt lJuvela et al.ll2012bt iPalmeirim et all 


20131) . none of these works to our knowledge has identified and 


analysed specifically non-fragmented filaments that would allow 
an assessment of the filament structure close to its initial condi¬ 
tion. 

Another key question is how the filaments fragment into 
dense cores that form stars. Analytic models predict the ini¬ 
tially hydrostatic, infinitely long filament to fragment through 
gravitational instability a t a wavelength of roughly four times 
the filament scale-height dlnutsuka & Mivam3ll992h . in a time- 
scale comgarablejojhesoui^-crossing time (order of a Myr, 
e.g., iFischera & Martini 1201 2ah . For finite-sized filaments, it 
has been hypothezised that the local collapse time-scale varies 
as a function of the position along the filament, and that the 
local c ollapse proceeds most rapidly a t the ends of the fil¬ 
ament (iBastienI 11983^ iPon et alT 1201 ih . While testing these 
predictions is an active goal of observational studies (e.g. , 
Jackson^t^ 20101: Schmalzl et al.l 2010t Busauet et al.ll2013l: 


Kainulainen et al.ll2013 : Takahashi et al.ll2013l) . it is of great in¬ 

terest to find and characterise the fragmentation in as simple as 
possible filaments to isolate different processes playing a role. 

The main technical difficulty in addressing the questions 
above is the general difficulty of tracing the distribution of 
molecular gas in the interstellar medium: all commonly-used 
tracers (dust emission, extinction, molecular line emission) 
have their specific problems and none of them can alone 
trace the wide range of column densities present in filaments. 
An additional difficulty is introduced by the small spatial 
scale of the filaments and their fragments. Filaments typically 
show widths that are on the order of 0.1 pc (e.g., iBarnardI 
119271: [Arzoumanian et al.ll201 Ih : fragments are similar in size. 
Consequently, single-dish observations can only study the struc¬ 
ture of filaments in adequate detail in the most nearby molecu¬ 
lar clouds (distances < 200 pc). As a result, it remains unknown 
what the physical structure of filaments is like prior to their frag¬ 
mentation and how exactly their fragmentation proceeds from 
that initial state to star formation. 

To make progress in characterising the early evolution and 
fragmentation of filaments, we target in this study the remark¬ 
ably long (~6.5 pc), highly rectilinear Musca molecular cloud. 
The cloud is located in the Musca-Cham aeleon cloud comple x 
(see Fig. [T]) at the distance of ~150 pc dKnude & HodflO*^ . 
The cloud has been identified to harbor only a very small num¬ 
ber of cold cores given its spati al extent (IVilas-Boas et al.iri994l: 
Lluvela et al.l2010l[201 lLl^i2bh . and it contains only one T Tauri 
star candidate, attesting to its low current star formation activ¬ 
ity (IVilas-Boas et aljri994t) . We have recently shown using CO 


line obse rvations that Musca is entirely subsonic and velocity- 
coherent dHacar et akllsubmittedh . Contrasting this apparent qui¬ 
escence, we have also shown that the probability distribution of 
column densities in the cloud shows a power-law-like excess 
of high-column d ensities, a feature common fo r star-forming 
molecular clouds dKainulainen et alJl2009ll20Tll) . Similarly, the 
cloud’s volume density distribution is more top-heavy than that 
of other nearby starless cloud s, suggesting that Musca may be 
close to active star formation dKainulainen et alJl2014ll . or pos¬ 
sibly in transition between non-star-forming and star-forming 
states. 

The above properties make Musca a perfect candidate to be 
in the earliest phases of fragmentation. In this paper, we present a 
detailed study of the column density structure of Musca. We use 
archival near-infrared (NIR) observations and 870 pm dust emis¬ 
sion observations to map column densities through the cloud and 
to characterise its internal structure. We then place the results 
in the context of hydrostatic equilibrium models, establish their 
suitability to describe the cloud, study its fragmentation, and fi¬ 
nally propose an evolutionary scenario for the cloud. 


2. Data and methods 


2.1. NIR data and dust extinction mapping 

The Musca c loud was observe d with the NIR wide-field camera 
NEWFIRM dAutrv et al.l2003h mounted on the Blanco 4-m tele¬ 
scope at the Cerro Tololo Inter-American Observatory between 
2011 April 11 and 13. We retrieved the raw data and associated 
calibration files from the NOAO public archive. Six pointings (5 
covering the entire cloud and 1 off the cloud) were obtained in 
the J, H and Kg broad-band filters. Each pointing was observed 
in dither mode. Table[T]gives a summary of the observation strat¬ 
egy. The instrument footprints are shown in Eig. [1] right panel; 
note that only four frames fall along the main part of Musca 
shown in the figure. 

The individual ra w images were p rocessed using an updated 
version of Alambic (IVandamel l2002h . a software suite devel- 
opped and optimized for the processing of large multi-CCD im¬ 
agers. Alambic includes standard processing procedures such as 
overscan and bias subtraction, flat-field correction, bad pixel 
masking, CCD-to-CCD gain harmonization, sky-subtraction, 
and de-stripping. Aperture and PSE photome try were extracted 
from the individua l images usin g SExtractor (iBertin & ArnoutsI 
1 19961) and PSEEx (iBertinI 120111) . The individual catalogs were 
then astrometric ally registere d and aligned on the same pho¬ 
tometric scale using Scamp (Bertin 20061). an d the final mo¬ 
saic produced using SWarp (iBertin et al.ll2002l) . The photomet¬ 
ric zero-points of the merged catalog were derived using the 
2MASS catalog. The final absolute astrometric accuracy is ex¬ 
pected to be better than O'.'! and the absolute photometric ac¬ 
curacy (rms of the individual images zeropoints) of the order of 
5%. The NEWEIRM data were complemented w ith J, H, and Kg 
band photometric data from the 2MASS archive dSkrutskie et al.l 
l2006h to cover a wider field towards Musca. The area covered by 
the NEWEIRM frames cover the area used in all analyses of this 
paper; the 2MASS data is merely used to complete the outskirts 
of a wide-field extinction map (Eig.[Tl right panel). 

We used the JHKs band data i n conjunction wi th the nicest 
color-excess mapping technique (iLombardil l2009t) to estimate 
the dust extinction throug h the Musc a cloud . Eor the details of 
the technique we refer to ILombardil (l2009l) . The technique is 
based on comparing NIR colors (H - K and J - H) and sur¬ 
face density of the stars that shine through the cloud to those 
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Fig. 1. Left: Dust extinction map of the Chamaeleon-Musca region, based on 2MASS data dKainulainen et al.ll2009ll . Right: Zoom- 
in to the Musca cloud, showing the instrument footprints of the observations employed in this work: NEWFIRM/NIR data (red 
rectangles), LABOCA/870 /rm cont inuum emission (blue dashed lines, approximate coverage), SHeFl/C'^O line emission (green 
crosses, from iHacar et al.llsubmittedi) . 


Table 1. Overview of NFWFIRM observations 


Band 

Number of dithered positions 

Number of coadds 

Single coadd exp. time 

J 

20 

2 

30 s 

H 

25 

16 

4 s 

Ks 

30 

12 

5 s 


from a field that is assumed to be free from extinction. We note 
that since the photometric depths of the NFWFIRM and 2MASS 
data are different, these parameters were determined separately 
for the two data sets. In implementing nicest, a relative extinc¬ 
tion law between ba nds has to be assume d. We chose the 
empirical extinction law of ICardelli et al.l (Il989h 


tk = 0.60 X th = 0.40 X Tj. (1) 

We emphasize that the technique estimates the dust extinction 
in NIR. We convert the estimates of NIR extinction to visual 
extinction by adopting t he ratio between op tical depths in NIR 
and visual wavelengths dCardelli et al.lll98^ 


Tv = 8.77 X tk. 


( 2 ) 


The visual extinctions are finally used to estimate the total hy¬ 
drogen (ato mic + molecular) col umn densities using the ratio 
mea sured by Savage et al.l dl977h andiBohlin et aP d l978h (see 
also lVuong et al.ll2003 : Giiver & Ozelll200'^ 


Ah = 2A(H2) -h A(H) = 1.9 x 10^' cm^^ x Ay mag^ 


(3) 


The final extinction map of Musca is presented in Fig. ^ with 
higher-detail blow-ups shown in Fig. [3] The map has the reso¬ 
lution of 45 arc sec and uncertainty at low column densities of 
cr(Av)=0.2 mag. We analyse the column and volume density 
structure of the cloud using this map in Sections [3. ll and lJSl 


2.2. LABOCA dust emission data 

In addition to extinction maps, we carried out in May 2010 
870 pm continuum observatio ns of the Musca cl oud using the 
FABOCA bolometric camera dSiringo et al.ll200^ at the APFX 
telescope. The observations consisted of fourteen individual 
raster-spiral maps combined to obtain a large-scale mosaic of 
approximately 2.2 by 0.4 deg^ covering the total extension of 
this cloud. Guided by o ur previous extinction maps of Musca 
dKainulainen et al.ll2009l) . the location, coverage, and overlap of 
these submaps were selected to optimise the sensitivity of the 
final mosaic along the main axis of the cloud. Skydips plus flux 
calibrations and pointing corrections in nearby sources were car¬ 
ried o ut every 1.5-2 hours. Using the BoA software dSchulleil 
1201 2h . the da ta reduction process fo llowed the standard method 
described by iBelloche et al.l d201 ih in order to recover the ex¬ 
tended emission in the cloud. Convolved to a final resolution of 
30 arcsec, the resulting map presents a typical rms of 17 mJy 
beam * and is shown in Fig]^ 

2.3. SHeFI and tine emission data 

We employ in the discussion of this paper *^CO (7 = 2- 1) and 
(7 = 2 - 1) line emission data of the Musca cloud from 
iHacar et^ dsubmittedh . to which we refer for the description of 
the observations and data reduction. The data were gathered with 
the SHeFI instrument at the FSO/APFX telescope in Chile. The 
spatial resolution of the data is 2875 and the velocity resolution 
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0.10 km or 0.17 km (the receiver backend was changed 
during the observations). The spectra typically have a main beam 
temperature rms of 0.14 K. Line emission was observed in 305 
positions along a curve that approximately traces the crest of the 
Musca cloud (see Fig. [TJ. The data were used to determine the 
central velocities, Vlsr, and velocity dispersions, cry, along the 
filament. They were also used to estimate the non-thermal con¬ 
tributions to the observed line-widths, ctnt and ctnt.coit (opacity- 
corrected), assuming a temperature of 10 K. For the details of the 
derivation of these values we refer to lHacar et alJ dsubmittedh . 

3. Results 

3.1. Column density structure of the Musca filament 

Figure |2] shows the column density map of the Musca cloud de¬ 
rived from the NIR dust extinction data. Figure |3] shows zoom- 
ins of the map. The highest column densities in the map are 
about 35 X 10^' cm~^. The mean column density above 2 x 10^' 
cm“^ is <A^h) = 4.5x10^^ cm“^. Figure|2]also shows the contours 
of the LABOCA 870 jum emission data. The emission contours 
show good morphological correspondence with the NIR-derived 
column density map, evidencing that no column density peaks 
remain undetected from the NIR-derived data, e.g., because of 
the saturation of extinction or resolution effects. In general, the 
relationship between the extinction- and emission-derived col¬ 
umn densities shows a high scatter, likely due to the spatial fil¬ 
tering of the dust emission data and their lower signal-to-noise 
compared to extinction data; the relationship is presented in 
AppendixiAl 

We estimated the total mass of the cloud from the NIR- 
derived column densities by summing up all column densities 
inside the contour^ of Ah = {1.5,3,4.5) x 10^' cm“^ (see 
Fig. [3] for the contours). This resulted in the total masses of 
M - {281,174,132) M©. We used the total masses to calcu¬ 
late equivalent line masses, M®, of the entire Musca cloud. The 
length of Musca for all thresholds is about I - 6.5 pc, resulting 
in My - {43,27,20) Mq pc"*. The inclination angle, i, of the 
cloud is not known, and for reference, we adopt here i - 0. The 
inclination affects the line masses through a factor of cos i, and 
hence, the values we give represent upper limits. Note that in ad¬ 
dition to the inclination uncertainty, the uncertainty of the line 
mass is significantly affected by distance uncertainty, which is at 
least some ~10%. 

The column density map (Fig. |2]i reveals numerous details 
within the large-scale Musca cloud. We divided the cloud into 
three regions, namely South, Center, and North, based on their 
morphology. The detailed column density maps of these regions 
are shown in Fig. [3] 

The column density map (and emission data) show that 
Musca contains several small-scale condensations. To have a 
simple, quantitative census of these structures, we identify “frag¬ 
ments” in the observed region. We emphasize that our purpose 
here is merely to have a first look at the small-scale substructures 
in the cloud; a simple fragment definition will suffice for this 
purpose. With this in mind, we first smooth the resolution of the 
emission data down to the resolution of the extinction data (45") 

* To describe how the chosen contour level affects the results, we 
give the masses and line masses corresponding to three different contour 
levels. 

^ We use the term “equivalent line mass”, because the cloud is partly 
fragmented (North and South regions) and does not resemble a well- 
defined filament throughout its length. In this case, “line mass” alone 
would be a somewhat misleading term. 


and identify fragment candidates as local emission maxima. In 
identifying the local maxima, we use contour steps of 0.03 mJy. 
If the fragment candidate is smaller than 45" in diameter, it is re¬ 
jected. As a second step, the fragment candidates are inspected 
against the local maxima of the extinction map. Local maxima 
are determined from the extinction data using contour intervals 
of 1.5x10^' cm"^. If the lowest contour that outlines a fragment 
candidate does not overlap with the lowest contour outlining an 
extinction maximum, the candidate is rejected. The pixel with 
maximum emission value is recorded as the location of the frag¬ 
ment. As a result of these steps, our fragments are essentially 
local emission maxima that have counterparts in the extinction 
map. The procedure results in identifying 21 fragments from the 
observed area. 

We estimated the masses associated with the fragments by 
summing up all column densities from the NIR-derived data 
within circular apertures. The radii of the apertures, listed in 
Table |2] were chosen by eye. The fragment masses span 0.3- 
7.1 Mq (Table |2 ]i. We note that all column densities within the 
apertures were included in the calculation, i.e., possible ’’back¬ 
ground” was not considered; it is not our goal to analyse the 
internal structure of the fragments or their relation to their sur¬ 
roundings. Finally, we roughly estimated the mean densities of 
the fragments from their area and mass, assuming spherical sym¬ 
metry (Table |2]l. 

In the following, we describe the specific properties of the 
South, Center, and North regions separately. The South region is 
clearly fragmented into five core-like structures, the most mas¬ 
sive of which contain several solar masses. We measured the sep¬ 
aration between the fragments as distances between neighbour¬ 
ing fragments; in case of the South (and Center) this is trivial, 
because the fragments are arranged into a line along the fila¬ 
ment crest. The mean separation in South is 0.78 + 0.6 pc. Note 
that the mean is strongly affected by the long distance between 
fragments #3 and 4; without the separation between those two 
fragments the mean peak separation is 0.50 + 0.2 pc. The total 
mass of the South region above Ah = {1.5,3,4.5) x 10^* cm"^ 
are M = {96,88,74) Mq. The region is about 2.7 pc long, result¬ 
ing in the equivalent line masses of = {36,33,27) Mq pc"*. 
Finally, we note that fragments #4 and #5 practically co mpose 
the object studied as “a filament” bv iJuvela et al.l (l2012bl) . It ap¬ 
pears to us that this region is clearly fragmented, and therefore, 
the axisymmetric cuts across the region measure in effect a com¬ 
bination of the profile of the fragments and that of the collapsing 
filament, not the profile of the filament itself. 

The Center region of the Musca cloud consists of a 1.6 pc 
long filament that shows less sub-structure than the South re¬ 
gion and has a relatively constant peak column density of Ah ~ 
10- 15 X10^* cm"^ along its crest. The line mass of this filament 
using the column density thresholds of Ah = {L5,3,4.5) x 10^* 
cm"^ results in M\ - {31,26,21) Mq pc"*. Thus, its line mass 
is practically the same as that of the entire Musca cloud on av¬ 
erage and that of the South and North (see below) regions. The 
FWHM of the filament is 0.07 pc, measured from its mean col¬ 
umn density profile (see Section l3^ . There are five fragments 
along the main filament and one about 0.5 pc off the filament. 
The fragments along the filament have very low contrast to their 
surroundings compared to those in the South (see the quan¬ 
tification of the column density contrast later in this Section). 
The mean separation of the five fragments along the crest is 
0.35 + 0.13 pc (fragment #9 was not included in the analysis). 
We consider physical models for the Center filament in Section 

im 
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Fig. 2. Column density map of Musca (colour scale), derived using near-infrared dust extinction mapping and data from NEWFIRM 
and 2MASS. The spatial resolution of the map is FWHM - 45" and the sensitivity cr(N}i) ^ 0.4 x 10^' cm“^. The white boxes 
indicate three regions that are shown in blow-ups in Fig. [3] The contours show the 870 /rm continuum emission contours, observed 
with the APEX/FABOCA instrument and smoothed down to the resolution of the extinction data. The contours are drawn at the 
intervals of 0.03 mJy/beam, which is the 3-crrms noise of the smoothed data. 


The North region of the Musca cloud is clearly fragmented, 
similarly to the South region. In the south-west corner of the re¬ 
gion, there is a complex area containing four fragments. This 
region is characterised by two velocity components, likely re - 
sulting from two overlapping structures (iHacar et ahllsubmittedh . 
In the middle of the North region, there is a curious structure that 


consists of two thin (Ah ~ 3 x 10^' cm“^) structures parallel to 
each other (a bubble-like structure), and of two peaks (fragments 
#19 in northwest and fragments #17 and 18 in southeast) that are 
elongated almost in perpendicular direction to those two struc¬ 
tures. Unfortunately, we did not cover this region with molec¬ 
ular line data, and thus the nature of the structure remains un- 
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Fig. 3. Blo'w-ups from the NIR-derived column density map of the Musca cloud (Fig. |2]l. The contours start at A^h = 1.5 x 10^* 
cm“^ and the step is 1.5 x 10^* cm“^. The lo'west contour is dra'wn 'with a dotted line. The dashed lines in the Center region indicate 
the 1.6 pc velocity-coherent filament 'whose radial structure is modelled in Section [T2l The numbered signs refer to the fragments 
identified 'within the cloud (see Section lTTI) . 


clear. There is a strong local maximum at the north-east corner 
of the region. This peak coincides 'with the loc ation of the only 
candi date young star (T Tauri star candidate, IVilas-Boas et alJ 
1 19941) in the cloud. The North region harbours ten fragments 
(Table|2]i. These fragments are not aligned along a line as clearly 
as fragments in Center and South. Ackno'wledging this caveat, 
■we calculated the separations along t'wo routes that connect the 
fragments (in the order of the fragment number): one includ¬ 
ing fragment #19 and excluding fragments #17 and 18, and the 
other excluding fragment #19 and including fragments #17 and 
18. This resulted in mean separations of 0.24 + 0.16 pc and 
0.25 + 0.17 pc, respectively. The total masses of the region above 
Ah = {1 -5,3,4.5) X 10^' cm^^ are M = {59,47,37) Mq. With the 
approximate length of 2.0 pc of the region, these result in the 
equivalent line masses of M° = (30,24,19) Mq pc“^. 

We emphasize that the North and South regions contain high- 
contrast density enhancements compared to the smoother, fil¬ 
amentary Center region; even though 'we identify ’’fragments” 
from the Center region, their column density contrast against 
their surrounding, more elongated gas component is very lo'w. 
We further quantify this morphological difference by measur¬ 
ing the standard deviation of column density (extinction) values 
along the approximate crest of the cloud (see Fig.[T]right panel). 
The histograms of the extinction values along the crest in differ¬ 
ent regions are sho'wn in Fig. H) The figure sho'ws that the dis¬ 
tributions of North and South regions are clearly different from 
that of the Center region: the North and South peak at around 
Ah ~ 6 X 10^' cm“^ and have strongly ske'wed distributions, 
■while the Center peaks at around Ah ~ 13 x 10^* cm“^ and has 
a narro'w distribution. The relative standard deviations of the ex¬ 
tinction values along the cloud crest are: 55% for South, 49% 
for North, and 15% for the Center. This sho'ws that the column 
density variation in the Center is clearly smaller, by a factor of 
~3, than in the North and South. 


3.2. Physical models for the Center filament 

The column density maps of the Center region (Figs. l2l|3]l re¬ 
vealed an about 1.6 pc long, relatively non-fragmented filament. 
In this section, 'we examine the radial structure of this filament 
and consider physical models for it. 

We estimated the mean column density profile of the Center 
filament from radial cuts across it. For this purpose, 'we first de¬ 
fined the filament crest by connecting the local maxima 'within 
the Ah = 8 x 10^' cm“^ contour 'with lines. The filament crest is 
sho'wn in Fig.|B]and the points defining it are given in Table lB.11 
We then sampled the radial profile 'with cuts perpendicular to the 
crest. The positions of the cuts along the crest 'were separated by 
4 pixels to ensure that the profiles at the crest position are not 
correlated. Consequently, the 1.6 pc long filament 'was sampled 
■with 20 profiles. 

Figure |5] sho'ws the resulting mean profiles, given separately 
for the 'west and east sides of the filament. Importantly, the pro¬ 
files sho'w that the filament is not symmetric, but the 'west side of 
it is steeper than the east side. The CO spectra sho'w t'wo velocity 
componen ts on the east side of the filament crest (cT. lHacar et al.l 
Isubmittedh . suggesting that another, faint structure overlaps 'with 
the main crest at the east side. This asymmetry is also seen at 
larger scales in the extinction map of the region (see Fig. [T]i. 
Therefore, 'we consider that the 'west side profile reflects the 
true filament profile more correctly. It is note'worthy that even if 
Musca seems to represent one of the simplest filamentary clouds 
■we kno'w, its radial profile is asymmetric 'with respect to the main 
axis 'when resolved in high detail. Similarly, as sho'wn in Section 
13.11 small scale variations are present even in the simple Musca 
filament. These phenomena are likely present in more complex 
filam entary clouds, suc h as L1495-B213 (e.g., ISchmalzl et al.l 
I 2 OIOI) and IC5146 ('e.g.. l7krzoumanian et aljl201 ll) . This advices 
caution 'when analysing profiles of complex filaments by averag¬ 
ing individual profiles and need to be addressed for the correct 
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Table 2. Physical properties of fragments in Musca 


# 

R.A. 

Dec. 

[10^‘ cm-2] 

T870,peak 

[Jy beam^*] 

R 

[pc] 

3^nir 

[Mo] 

(n) 

[lO'* cm-3] 

Region 

1 

12:21:26 

-72:21:06 

23 

0.15 

0.07 

2.5 

4 

South 

2 

12:21:54 

-72:17:02 

19 

0.11 

0.07 

1.7 

3 

South 

3“ 

12:22:19 

-72:10:24 

17 

0.10 

0.06 

1.0 

2 

South 

4 

12:24:32 

-71:50:23 

27 

0.22 

0.08 

3.5 

3 

South 

5 

12:25:37 

-71:42:14 

34 

0.18 

0.11 

7.1 

3 

South 

6 

12:26:46 

-71:30:00 

15 

0.10 

0.07 

1.7 

3 

Center 

7 

12:26:59 

-71:27:41 

17 

0.10 

0.06 

1.2 

3 

Center 

8* 

12:27:45 

-71:22:07 

15 

0.10 

0.05 

1.3 

3 

Center 

9" 

12:29:46 

-71:22:42 

6 

0.05 

0.04 

0.5 

1 

Center 

10" 

12:28:30 

-71:18:13 

17 

0.12 

0.04 

1.4 

4 

Center 

11“ 

12:29:26 

-71:10:41 

20 

0.17 

0.05 

1.4 

3 

Center 

12 

12:30:48 

-71:10:08 

15 

0.14 

0.05 

0.9 

4 

North 

13 

12:31:03 

-71:04:32 

19 

0.10 

0.04 

1.1 

6 

North 

14 

12:31:05 

-71:01:54 

19 

0.09 

0.04 

0.9 

5 

North 

15 

12:31:53 

-71:02:08 

17 

0.07 

0.07 

1.8 

3 

North 

16 

12:32:13 

-70:57:49 

13 

0.06 

0.04 

0.6 

4 

North 

17 

12:34:09 

-70:57:40 

7 

0.06 

0.04 

0.3 

2 

North 

18 

12:34:24 

-70:56:11 

8 

0.08 

0.04 

0.3 

2 

North 

19 

12:32:41 

-70:44:21 

15 

0.19 

0.06 

1.0 

2 

North 

20 

12:35:00 

-70:43:54 

11 

0.12 

0.05 

0.6 

3 

North 

21 

12:35:17 

-70:40:10 

30 

0.62 

0.05 

1.1 

5 

North 


Notes. The peak value of dust emission from the fragment is more than 45"ofFset from the associated extinction peak. The structure is 
considered as a fragment, because the contours outlining the local maximum in emission overlap with the contours outlining a local maximum in 
extinction. The mass estimate is based on extinction data; the aperture for the mass calculation was centered on the nearest extinction maximum. 

The contours defining the local maximum in emission map overlap with two local extinction maxima. The mass estimate is based on extinction 
data; the aperture for the mass calculation was centered on the nearest extinction maximum. The fragment is offset by about 0.5 pc from the 
filament crest. Consequently, we do not include the fragment into the fragmentation analysis. 


physical interpretation of the derived parameters. In the follow¬ 
ing, we quantify the radial structure of the Center filament in 
Musca by fitting different physically motivated models to its ob¬ 
served profile. 


3.2.1. Isolated, infinite filament in hydrostatic equilibrium 

We first modelled the 1.6 pc long filament in the Center region 
in the context of an infinite, isolated cylin der in equilibriu m be¬ 
tween thermal and gravitational pressures (IOstrikeilll964ll . Such 
an equilibrium configuration has the line mass 




2kT 

finiiiG 


2 ^ 'T 

—^ 16 Mq pc“' X (-), 


(4) 


where Cs is the sound speed. This indeed is close to the ob¬ 
served line mass of the Center filament (Mi = 21-31 M©). We 
note that it is possible that non-thermal motions provide addi¬ 
tional support for a filament, making its critical line mass higher. 
An estimate of this “turbulent support” can be calculated from 
Eq.m by replacing the sound speed with an effective sound 
speed, Cs,eff, that results from the total line width in the cloud. 
The non-thermal line-width of C'^O) is about 0.7cs, yielding 
Cs.eff = Vl -H O.T^Cs = 1.22cs. The critical line mass is affected 
by the square of this, yielding 24 M© pc“*. 

The volume densit y profile of the infinite, hydrostatic fila¬ 
ment in equilibrium is (IOstrikeilll964li 


p(r) = 


_Pc_ 

(l+(r/(V8ro))2)/’/2’ 


(5) 


where ro is the scale radius as defined bv lOstrik^ (1 19641) . pc the 
central density, and the exponent p=4. Note that if the exponent 


is not 4, the radial profile does not describe an equilibrium solu¬ 
tion. Using the generic, Plummer-like profile is however useful, 
because it allows for quantification of ro, pc, and p through a 
fit to an observed profile, and from therein, examination of the 
validity of the model. The scale radius is coupled to the central 
density and temperature, T, through 


pm-aAnGr^ 

The column d ensity profile correspondin g to the volume density 
profile is (e.g.. I Arzoumanian et al.ll20Tll) 


N{r) = Ap 


_ Pc Vs^o _ 

(l+(r/(V8ro)f/^'^^^" 


(7) 


In the equation, Ap = (i+^pp ’ i as the inclination 

angle of the filament, which we assume to be zero. 

Figure |5] shows the best fit of Eq. |7] to the observed mean 
profile, achieved by treating pc, ro, and p as free parameters, 
fixing the fitting range to r = [0,0.6] pc, and then minimising 
the chi-square between the data points and the model. The in¬ 
verse variances of the 20 individual profile cuts were used as the 
weights of the data points in the fit. The best-fit model for the 
west-side profile has the exponent p = 2.6+11%, the scale radius 
ro = 0.016 ± 33% pc, and the central density = 54000 + 50% 
cm 2 (Table |3]l. The quoted errors are the one-sigma uncertain¬ 
ties of the fit. One should note that fitting Eq. |7] in this man¬ 
ner suffers from degenera cy between pc, ro, and p (lJuvela et al.l 
12012at ISmith et al.n2014l) . This degeneracy likely causes the 
large uncertainties in ro and pc. The slope p suffers less from 
the degen eracy and in gener al can be constrained with lower un¬ 
certainty dSmith et al.ll2014h . 
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Fig. 5. Mean radial column density profiles along the Center filament of the Musca cloud. The left and right panels show the mean 
profile on the west and east sides of the filament axis, respectively. The blue line shows the fit of Plummer-like profile, and red line 
shows the the fit of a pressure-confined hydrostatic cylinder. The colored arrows indicate the radial ranges used in fitting the profiles. 
The error bars show the standard deviation of 20 individual profiles from which the mean profile is constructed. The vertical dashed 
line shows the FWHM size of the extinction mapping beam. 

Table 3. Parameters of the best fit radial density distribution models. 


filament model 

AC 

[pc] 

Rf 

[pc] 

P 

ro 

[pc] 

tic 

[10^ cm-3] 

/cyl 

Pext/fc 

[lO'* Kcm-3] 

other 

Plummefi’, west 

<0.6 

- 

2.6±11% 

0.016±33% 

5 

- 

- 

- 

Plummer* , east 

<0.6 

- 

1.8±5% 

0.012±38% 

5 

- 

- 

- 

P-confined'^ 


0.12±7% 

- 

- 

10 

0.66±12% 

12±37% 

- 

P-confined, with 

<0.6 

0.10±12% 

- 

- 

10 

0.64±18% 

13±52% 

Re=0.6 pc, No= 


envelope'^ 2.4x10^' cm ^ 


Notes. The fitting range. Eq.[7l Eq.[8] Eqs.[8]and|^ 



N [10^' cm^G 

Fig. 4. Frequency distribution of column density values along 
the approximate crest of the Musca cloud for South (red). Center 
(black), and North (blue) regions. The dashed lines indicate the 
mean column densities. The relative standard deviation of the 
column density values in each region is shown in the panel. 


3.2.2. Isothermal filament in equilibrium, confined by external 
pressure 


We also modelled the filament in the Center region as a pressure- 
confined (an d truncated), isothermal f ilament in hydrostatic 
equilibrium (iFischera & Martini 12012ah . The column density 
profile of this model is 


^ (A'^h) ' 


4nG 1 - /cyi 

1 “ /cyl 


1 /cyl + ^ /cyl 


7/cyl(l-/cyl)(l-x2) 


1 /cyl 


l-/cyl(l-^2) 


X arctan 


/cyl(l - X^) 




(8) 


In the equation the dimensionless coordinate x is the radius, r, 
normalized by the outer radius of the filament, Rf, i.e., x = r/R{, 
and /cyl = Mi/Mi,crit is the ratio of the filament’s line mass to the 
critical line mass, (see Eq.|4]for Mi crit)- The relative shape of the 
column density profile depends only on the parameter /cyi and 
on the normalisation of the radial profile, i.e., on Rf. The tem¬ 
perature and line mass affect the profile shape implicitly through 
/cyl. In addition, the normalisation of the profile is affected by 
the confining external pressure pext- 
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The applicability of this model can be immediately examined 
through /cyi, which in the context of this model cannot be larger 
than unity. In Section lTTl we estimated the line mass of the fila¬ 
ment to be Ml = {31,26,21) Mq when using the column density 
thresholds of A^h = {1-5,3,4.5) x 10^' cm“^ to compute the line 
mass. The critical line mass without the contribution from non- 
thermal motions (16 Mq pc“') is slightly below this range, and 
the critical line mass with the non-thermal contribution (24 M© 
pc ') overlaps with it. One should note that the observed line 
masses are upper limits because of the assumption of zero incli¬ 
nation. Distance uncertainty is also at least some 10%. We also 
do not know exactly the temperature structure of the filament, 
which also affects Mi crit- Finally, in the context of this model, 
the filament is pressure-confined by the surrounding gas; this 
surrounding gas possibly contributes to the observed Mi value 
(we address this further below). We conclude that the observed 
line masses are such that considering the pressure-confined fila¬ 
ment model is reasonable. 

We fitted Eq. [8] to the observed data points at radii below 
0.6 pc by letting /cyi and pext change as free parameters, again 
using the inverse variances of the 20 profiles to weight the data 
points. This did not yield an acceptable fit. However, allowing 
also the fitting radius, i.e., the radius of the filament, to change 
as a free parameter and ignoring the profile outside it resulted 
in a reasonable fit (see Fig. |5]). The best fit has the parameters 
/cyi = 0.66 + 12%, and pext = 12 + 37% x 10"^ K cm“^, and Rf - 
0.12 + 7% pc where the uncertainties result from the chi-square 
fit (Table O. The external pressure is consistent with pressures 
needed to confine isolated gl obules modelled as Bonnor-Ebert 
spheres (~2-20xl0^ K cm~^. l^ndori et ^l2005h . and a factor 
of a couple hig her than inferred for parse c-scale filaments (~1- 
5x10^ K cm- MFischera & Martinll2012all^ . 

One has to recognise that the above model does not include 
any contribution to the observed column densities from the sur¬ 
rounding gas that is r equired to provide the confining pressure 
(see the discussion in iRecchi et al.ll20T4l) . Thus, the use of Eq. 
[8] by itself may not be physically entirely meaningful. In prac¬ 
tice, the confining pressure could arise from the thermal pressure 
in a two-phase medium, in which the cooler, denser component 
is surrounded (and confined) by the less dense, warm compo¬ 
nent. Accretion of gas in to the filament could also possibly pro¬ 
vide confining pressure (lHeitschll2013alf5l . Turbulent ram pres¬ 
sure could also provide confining pressure, although to provide 
confinement, it should be relatively isotropic, which may be un¬ 
realistic. Regardless of the nature of the pressure, a surrounding 
gas component is needed, and it can affect the observed column 
density structure. 

To explore the effect of the surrounding gas to the parame¬ 
ters of this model, we made an experiment in which the pressure- 
confined filament is embedded inside a cy lindric envelope (e.g., 
iHernandez et alJl2012HRecchi et alJl2Q14h . This cylindric enve¬ 
lope has a power-law radial density distribution, Pe(r) = Cr^^, 
where C is a proportionality constant. The column density pro¬ 
file of the envelope res ults then from integ rating Pe(r) along the 
line of sight (see, e.g., iRecchi et alJl2014ll . yielding 

arctan ■\J(Relr)^ - 1 

- arctan -\/{R{lr)^ - 1, r < Rf (9) 

^arctan ■\J(Relr)^ - 1, Rf < r < R^, 

where Re and Rf are the radius of the envelope and filament, 
respectively. At the limit r —> 0 the equation has the value 
C(-Rjr* + Rf *), which gives the column density of the envelope 
at r = 0, i.e., Nq. The total column density of the model is the 



r [pc] 


Fig. 6. Mean radial column density profiles along the Center fil¬ 
ament. The red line shows the fit of a pr essure-confined equi¬ 
librium model (iFischera & Martini 12012ah . embedded inside a 
cylindric envelope. The black dashed line shows the column den¬ 
sity profile of the filament without the surrounding envelope. The 
dashed vertical line shows the FWHM size of extinction map¬ 
ping beam. 

sum of the envelope and the embedded filament (Eq. [^i. This 
model was fitted to the observed profile with /^yi, pext. Rf, and 
C as free parameters, assuming that the envelope goes to zero at 
0.6 pc, i.e.. Re = 0.6 pc (Fig. |6]l. We stress here that our goal 
is to probe how the filament parameters are affected by an en¬ 
velope, not to find the best possible envelope model, nor even a 
physically meaningful one. It is evident that the envelope model 
we use is rather arbitrary, and other envelopes may well produce 
equally good fits. 

We find that the observed profile is best fit when a filament 
is surrounded by an envelope that has the central column den¬ 
sity No = 2.4 X 10^' cm“^. The parameters of the embedded 
filament, i.e., /cyi and pext (see TableO are not significantly dif¬ 
ferent from those of the filament without the envelope. In con¬ 
clusion, a pressure-confined hydrostatic filament, embedded in a 
low-density envelope, seems to provide a reasonable explanation 
for the entire observed column density profile. 

3.3. Relationship between the column density and velocity 
structure in Musca 

We briefly compare in this section the column density struc- 
ture of the Musca clou d to the large-scale kinematic data from 
[Hacar et^dsubmittedl) . While Hacar et al. discusses the proper¬ 
ties of the global velocity field in Musca, we put here those data 
in the context of the column density structure and fragments of 
the cloud. 

Figure |7] shows the NIR-derived column densities, '^CO and 
C'^O central velocities (Visr), non-thermal velocity dispersions 
(c^nt), and non-thermal opacity-corrected velocity dispersions 
(cNT.corr) for the ~300 line-of-sights observed along the Musca 
cloud (see Fig.[T]). Note there are some double-peaked spectra at 
the edges of the cloud; Fig. |7] shows the Visr and cnt values for 
both these comp onents, obtained via fit s of separate Gaussian 
components (see iHacar et alJ[submitted for details). The loca¬ 
tions of the fragments that are closer than 100 arcsec of the cloud 
crest are al so shown in the figure. As found out by [Hacar et al.l 
dsubmittedh . Musca is velocity-coherent and subsonic at densi- 
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N,(r) = - 
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ties probed by the C'^O data (sub-/transsonic at densities probed 
by the '^CO data) practically throughout its length. The veloc¬ 
ity gradient along the South region is about dVisr/dr « 0.5 km 
s ' pc '. The South r egion also shows local velocity gradients, 
up to ~3 km s ' (cf.. iHacar et akllsubmittedh . that clearly coin¬ 
cide with the column density fragments. The North and Center 
regions show a velocity gradient of dV^r/dr » 0.15 km s“' pc“', 
with the only local gradients coinciding with the densest frag¬ 
ment in the North region (#21) and the complex clump of multi¬ 
ple maxima (close to fragment #13). In particular, the two most 
pronounced fragments in the Center region (#6 and 7 in Table 
13 do not show gradients in C'^O central velocities. The scatter 
of the line velocities in The North and South is clearly larger 
than in the Center (dominantly due to the large-scale gradient): 
(r(yLSR,C'*^0),o„th = 2.03 km S-', fr(yLSR, c''0)ee„tre = 0.67 
km s ', cr(yLSR, C'^0)north = 1-57 km s '. The North and 
South also show slightly larger scatter of non-thermal velocity 
dispersions than the Center; cr((TNT, C'®0)south = 0.3 km s“', 
Cr(crNT, C"'0)centre = 0.2 km S"', Cr(crNT, C'*0)north = 0.24 km 
s"'. Finally, some oscillatory motions are seen at the locations of 
the two most prominent fragments in the South region (#4 and 
5). Their presence sugg ests that fragmenting m odes indeed are 
present in the cloud (cf.. [Hacar & Tafan3l201 Ih . 

Summarising, the above properties demonstrate the fact that 
even if Musca as a whole is subsonic and velocity-coherent, the 
opposite ends of it show slightly more active internal dynam¬ 
ics than the extremely quiescent Center filament. This activity is 
likely related to the fragments present in the ends of the cloud. 


4. Discussion 

4.1. Musca: caught in the middle of gravitational 
fragmentation 

What do our observations tell about the nature of Musca and its 
fragmentation? We have demonstrated that the physical struc¬ 
ture of the extremely quiescent Center filament of Musca is in 
a good agreement with the predictions for the line mass, ra¬ 
dial density distribution, and velocity structure of the pressure- 
confined hydrostatic equilibrium model. The cloud contains 21 
fragments (according to our definition of fragments), with the 
higher-contrast ones located in the South and North regions. 
Specifically, the North and South have a factor of three higher 
column density variance compared to the Center region. Some 
of the fragments contain several solar masses, and thus have the 
potential to be pre-stellar cores. Importantly, the equivalent line 
masses of the North and South are equal with the line mass of the 
Center filament, giving rise to a hypothesis that they may have 
fragmented from an initial condition that resembles the relatively 
non-structured Center filament. 

To further examine this hypothesis, we consider the pre¬ 
diction of the gra vitational fragmentation m odel (the pressure- 
confined filament, iFischera & Martinll2012ah for the fragmenta¬ 
tion length and time-scale of a cylinder that has the properties of 
the Center filament. In this model, filaments with line masses 
less than the critical line mass can fragment when perturbed. 
Filaments with line masses larger than critical (that are not equi¬ 
librium solutions to begin with) should collapse into a spindle 
rapidly. The wavelength of the faste st growing perturbations i s 
between Tm = 5 - 6.3 x FWHMcyi (IFischera & Martinll2012ah . 
The lower limit corresponds to the case fiyi 1 and the upper 
limit to /cyi —> 0. We measured the FWHM of the Center fila¬ 
ment to be 0.07 pc, which results in the fragmentation scale of 
/Ini ai 0.35 - 0.44 pc. The observed mean peak separation of the 


fragments in all regions is 0.41 + 0.37 pc and the median is 0.34 
pc (mean is 0.32 + 0.18 pc and median 0.33 pc without the long 
separation between fragments #3 and 4), overlapping with the 
predicted interval. The results are practically the same regard¬ 
less of how we compute the separations in the North region (see 
Section [3TT) . We also note that fragment #9 is not included in 
this analysis because of its location off the filament axis. In con¬ 
clusion, it appears that also from the fragmentation scale point- 
of-view, the North and South regions may have originated from 
an initial condition that resembles the Center filament. 

The fragmentation time-scale of the infinite, 
pressure-confined cylinde r depends on fiyi, T, and pext 
(IFischera & Marti^l2012al) . and is between 


2.95(1 -/cyl)TKp ,/cyl^l 

4.08tkp ,/cyi —> 0, 


( 10 ) 


where tkp = sJkTKAngmuGpext) ~ 0.2 Myr. This gives the 
range ~ 0.2 - 0.8 Myr when using the parameters of the best¬ 
fitting pressure-confined model (Table [3]). We note the caveat 
that the above predictions for the fragmentation scale and time- 
scale are for infinite cylinders; it may be that the predictions dif¬ 
fer in finite geometry. 

Why has the fragmentation proceeded strongly at the ends 
of the Musca cloud, but not at its center? One possible mech¬ 
anism in the context of gravitational fragmentation is provided 
by finite-size effe cts in cylindri c ge ometry (as opposed to th e 
infinite models of lOstrikenI 19641 and iFischera & Martinll2012^ . 
In finite filaments, the acceleration of gas depends on the rel- 


(lBastienll983l; Pon et al.l201 1 

, 2012llClarke & Whitworthl2015l 

see alsolBurkert & Hartmann 

2004 for a studv in a sheet geom- 


etry). These higher accelerations may translate into shorter local 
collapse time-scales, and therefore, into format ion of fragments 
at the ends of the fil ament prior to it s centre dPon et al.ll2Ql iL 
I 2 OI 2 I) . For example, iPon et al] (1201 iL see their Figs. 5 and 6) 
suggests that the local collapse might occur a factor of two-to- 
three faster at the end s of the filament than at its center. Similarly, 
IClarke & Whitwo^ (120151) found in their numerical simula¬ 
tions that the first fragments form at the ends of the filament, 
regardless of its aspect ratio. While ther e are some quantita¬ 
tive di ffere nces between the predi ctions of IClarke & WhitworthI 
(I 2 OI 5 I) and iPon et alJ(l201 lLl2012l) . both models predict that long 
(but finite-sized) filaments fragment starting from their ends. 
Consequently, they lay out the basis for explaining the fragmen¬ 
tation at the ends of Musca as a result of the end-dominated col¬ 
lapse caused by higher accelerations, and therefore, shorter local 

collapse time-scales. _ _ _ 

In the models of iPon et al] (1201 iL 1201 2l) and 
IClarke & WhitworthI (1201 5l) . high-aspect ratio filaments 
collapse globally and the fragments formed at the ends of 
the cloud fall into its centre; this can be used to gain further 
insight into the evolution of Musca under this framework. 
The time-scales for the global collapse (time-scale for the end 
fr agments to reach the cloud centre) are (following the notation 
of IClarke & Whitworthll2015l) 


fda = (0.49 H- 0.26 A) (Gpo)■‘^^ 


and 


^nnn — 


0.44 A (Gpo)“‘^^ 
0.98 A'/2 (Gpo)-'^2 


A < 5 
A > 5, 


( 11 ) 


( 12 ) 


where G is the gravitational constant and po the central density. 
The aspect ratio of Musca depends on how we define its width. If 
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Fig. 7. Column density and velocity structure of Musca, measured approximately along the crest of the cloud (see Fig.[T]). The top 
panel shows the dust extinction from NIR data. The second panel shows the central velocities, Visr- The third and bottom panels show 
the non-thermal velocity dispersions, cnt (non-corrected) and (XNT.corr (opacity-corrected), normalised by the sound speed Cs. The 
arrows in the bottom panel indicate upper limits. The red and blue points show the data for '^CO (2-1) and C'*0 (2-1), respectively. 
The locations of the fragments that are closer than 100 arcsec of the cloud crest are indicated with thin dashed lines. The thick 
dashed vertical lines separate from each others the South, Center, and North regions. The thick dashed horizontal line indicates the 
sonic velocity dispersion and the dotted horizontal line the two times the sonic velocity dispersio n. The lower right corne r of the 
middle panel indicates with a line a velocity gradient of 3 km s '. The figure is a reproduction from iHacar et al.l (Isubmittedh . 


we adopt the ~0.6 pc extent of the Plummer profile (see Fig.|5l), 
the aspect ratio 6.5 pc / 0.6 pc 5 follows; if we adopt the 0.1 pc 
radius of the pressure-confined filament, the aspect ratio is 6.5 
pc / 0.2 pc 32. This interval of A = 5 - 33 is very conserva¬ 
tive, as it represents the extreme choices for the filament width. 
For the central density, we use 10^ cm“^. The global collapse 
time-scale for the chosen range of aspect ratios is then between 


~0.5-3 Myr, depending on the model. We showed earlier that 
the fragmentation time-scale of an infinite pressure-confined fil¬ 
ament is between 0.2-0.8 Myr. Even if these estimates are highly 
uncertain, they suggest that the fragmentation time-scale may be 
significantly shorter than the global collapse time-scale. This is 
in agreement with the fact that we do not observe Musca to be 
globally collapsing, but rather forming multiple fragments along 
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it. This, in turn, suggests that the age of Musca is rather closer to 
its fragmentation time-scale than its global collapse time-scale, 

i.e., only some tenths of a Myr. Alternatively, additional pro¬ 
cesses must be supporting Musca against global collapse. 

Instead of the gravitational focusing, another possibility to 
explain the differences in fragmentation along Musca may be 
differences in the initial physical properties along it. In principle, 
the f ragmentation time-scale an ti-correlates with the line mass 
(e.g., lFischera & Martinil2012ah : however as shown earlier, the 
equivalent line masses of the North and South are the same as 
the line mass of the Center, stating that they should collapse in 
the same time-scale. It could also be hypothesised that difference 
in the dynamics of the regions plays a role: the North and South 
regions show significantly higher global velocity gradients than 
the Center region, which might affect fragmentation in them. 
Whether such an initial velocity gradient could promote frag¬ 
mentation remains to be explored via theoretical works. Finally, 
the time-scale also correbtes w ith the size of the initial density 
perturbations dPon et alJl201 ih . It cannot be ruled out that the 
North and South regions would have had higher density variance 
along them compared to the Center, inducing faster fragmenta¬ 
tion. 

Ultimately, one would like to compare the time-scales de¬ 
rived above with other age estimators in the cloud. However, the 
only direct sign-post of age in Musca is the single T Tauri star 
candidate in the North region (fragment #21). Its presence in¬ 
dicates that fragmentation in the North region (and supposedly 
in South) must have taken place on the order of ~0.5 Myr ago 
(Class 0-1-1 protostellar lifetime. lOunham et al.l 120141) , with an 
admittedly large uncertainty. We showed earlier that the time- 
scale of gravitational fragmentation in Musca is 0.2-0.8 Myr, 
coinciding with this estimate. Further studies targeting, e.g., the 
chemical ages at the different locations along the filament would 
provide additional, independent measures of age. 

Altogether, the morphology of Musca, its fragmentation, ra¬ 
dial structure, and velocity characteristics suggest that gravita¬ 
tional fragmentation is taking place in situ in Musca, mould¬ 
ing it towards star formation. This is in a qualitative agreement 
with our ea rlier results regard i ng the cloud’s column density 
structure: in iKainulainen et alJ (l2009tl we analysed the proba¬ 
bility density function of column densities (A-PDF) in Musca, 
showing that it exhibits (at least a start of) a power-law-like tail 
at high-column de nsities, similarly w i th act ively star-forming 
clouds. Further, in IKainulainen et akl (12014|) we showed with 
the help of volume density modelling that Musca currently con¬ 
tains only a very small amount of gas (~6 Mq) at densities 
high-enough for star formation. From the point-of-view of its 
column and volume density structure, the cloud must become 
more condensed locally to form stars; this could happen via the 
on-going fragmentation and gravitational collapse. Musca rep¬ 
resents a perfect target for studying this: the simplicity of the 
cloud is ideal for isolating the different processes driving frag¬ 
mentation. 

5. Conclusions 

We analyzed in this work the column density structure and frag¬ 
mentation of the 6.5 pc long, rectilinear Musca molecular cloud 
that shows very little signs of on-going star formation. We ex¬ 
amined the structure of the cloud with the help of near-infrared 
dust extinction data and 870 jum dust emission data. We also ex¬ 
amined the large-scale kinematic st ructure of Musca using 
and C'^O line emission data (from lHacar et alJisubmitted) . Our 
conclusions are as follows. 


1. The column density data show that the Musca cloud is frag¬ 
mented into 21 fragments that have masses between ~0.3- 
7 Mq. The Center region of Musca harbours a remarkably 
well-defined, 1.6 pc long filament that contains only low- 
contrast fragments. The filament shows only 15% relative 
column density variations along its crest, in stark contrast to 
~50% variations in the more pronouncedly fragmented ends 
of the cloud. We emphasize that the relative column density 
variations provide a physically motivated, simple measure of 
how fragmented a filament is; a low degree of fragmentation 
is necessary for characterising the initial condition of frag¬ 
menting filaments. 

2. The line mass of the Center filament is 21-31 Mq pc ', i.e., 
close to the critical value of a cylinder in hydrostatic equi¬ 
librium. The equivalent line masses (total mass divided by 
length) of the North and South regions are practically the 
same. 

3. The radial density structure of the Center filament is well de¬ 
scribed with either i) a Plummer profile that has a power-law 
exponent of 2.6, scale radius of 0.016 pc (FWHM of 0.07 
pc), and central density of 5 xlO"^ cm“^, or ii) a hydrostatic 
filament that is confined by external pressure imposed by a 
surrounding envelope. In this case the ratio of the line mass 
to the critical line mass is 0.66, central density is 10^ cm^^, 
and external pressure is ~10^ K cm“^. We note that even 
this simple and well-resolved filament shows an asymmet¬ 
ric column density profile with respect to the main filament 
axis. This advices caution when analysing profiles of fila¬ 
ments that are more complicated than Musca. 

4. The fragmentation length of the cloud, ~0.4 pc, agrees with 
the prediction of the gravitational fragmentation of an in¬ 
finite, pressure-confined cylinder. The fragmentation time- 
scale predicted by the model is 0.2-0.8 Myr and the global 
collapse time of the cloud 0.5-3 Myr. 

5. The ends of Musca show somewhat more variation in their 
C'^O velocity structure than the extremely quiescent Center 
region, both in the central line velocities and the scatter of 
non-thermal line widths. The central line velocities peak at 
the locations of the fragments, suggesting motions associ¬ 
ated with the fragments. 

6. The above characteristics give rise to a hypothesis for the 
evolutionary scenario of Musca: an initially 6.5 pc long, 
close-to hydrostatic filament is in the middle of gravitational 
fragmentation. Fragmentation started a few tenths of a Myr 
ago, and it has so far taken place only at the ends of the 
cloud, possibly because of the shorter local collapse time- 
scales compared to the cloud centre (“gravitational focus¬ 
ing”). The relatively non-fragmented Center region may still 
represent the close-to hydrostatic initial condition of the fila¬ 
ment, regardless of whether the fragmentation at the ends of 
the cloud is driven by gravitational focusing or some other 
process. 
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Appendix A: Relationship between 

extinction-derived and emission-derived column 
densities 

We measured the column densities through Musca using two 
techniques: NIR dust extinction mapping and 870 fim continuum 
emission. We briefly compare here the column densities result¬ 
ing from them. The extinction was used to estimate the column 
densities as explained in Section 12.11 To estimate the column 
densities from emission, we used the formula: 

RFy 

A^laboca = ^(A.l) 

where Fy is the flux density, /? = 100 is the gas-to-dust ratio. By 
is the Planck function in temperature Tdust = 10 K, p = 2.8 is 
the mean molecular weight, k = 1.85 cm^ g ' the dust opacity 
(lOssenkopf & Henning|ll99^ . and Q the beam solid angle. 

Figure |AT| show the relationship between the column densi¬ 
ties from the two techniques. The relationship has a large scat¬ 
ter, likely resulting from the spatial filtering of the LABOCA 
data. The relationship indicates that no significant column den¬ 
sity peaks remain undetected by the NIR data; there is no signif¬ 
icant saturation of the NIR-derived column densities compared 
to LABOCA-derived ones. We illustrate this further by high¬ 
lighting in Fig. lA. li the relationship in the area of the most mas¬ 
sive fragment (#5) of the cloud. For this particular fragment, the 
NIR-derived column densities are systematically higher than the 
LABOCA-derived ones. This could be, e.g., due to misestima- 
tion of the temperature used in Ea. IA.il 


Appendix B: Definition of the Center filament crest 

To measure the radial profile of the Center filament, we defined 
its approximate crest with lines connecting nine points. The co¬ 
ordinates of these points are given in Table IBT] This crest was 
sampled with 20 equally-spaced cuts perpendicular to it. 
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Fig.A.l. Relationship between emission-derived 
(A(LABOCA)) and extinction-derived (A(NIR)) column 
densities. The red symbols show the data in the area of the most 
massive fragment in the cloud (#5). The dashed line shows the 
one-to-one relationship. 
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